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Bandwidth Efficient and Rate-Matched 
Low-Density Parity-Check Coded Modulation 

Georg Bocherer, Member, IEEE, Patrick Schulte, and Fabian Steiner 


Abstract —A new coded modulation scheme is proposed. At 
the transmitter, the concatenation of a distribution matcher 
and a systematic binary encoder performs probabilistic signal 
shaping and channel coding. At the receiver, the ontpnt of 
a bitwise demapper is fed to a binary decoder. No iterative 
demapping is performed. Rate adaption is achieved by adjusting 
the input distribution and the transmission power. The scheme 
is applied to bipolar amplitude shift keying (ASK) constellations 
with eqnidistant signal points and it is directly applicable to 
two-dimensional quadrature amplitude modulation (QAM). The 
scheme is implemented by using the DVB-S2 low-density parity- 
check (LDPC) codes. At a frame error rate of 10“^, the new 
scheme operates within less than 1 dB of the AWGN capacity 
|log 2 (l + SNR) at any spectral efficiency between 1 and 5 
bits/s/Hz by nsing only 5 modes, l.e., 4-ASK with code rate 2/3, 
8-ASK with 3/4, 16-ASK and 32-ASK with 5/6 and 64-ASK with 
9/10. 


I. Introduction 

Reliable communication over the additive white Gaussian 
noise (AWGN) channel is possible if the transmission rate per 
real dimension does not exceed the capacity-power function 

C(P)= ^l0g2(l + P/l) (1) 

where P is the transmission power and P/1 is the signal- 
to-noise ratio (SNR). An important consequence of 0 is 
that to achieve power and bandwidth efficient transmission 
over a substantial range of SNR, a communication system 
must adapt its transmission rate to the SNR. Contemporary 
standards such as DVB-S2 0 support input constellations of 
different sizes and forward error correction (FEC) at various 
code rates. The combination of a constellation size with a code 
rate forms a transmission mode. By choosing the appropriate 
mode, the system adapts to the SNR. In Fig. we display the 
operating points where a frame error rate (FER) less or equal 
to 10“^ can be achieved by bipolar amplitude-shift keying 
(ASK) constellations and the DVB-S2 low-density parity- 
check (LDPC) codes. Observe that the SNR gap between the 
operating points and the capacity-power function varies with 
the SNR. At the upper comer points, the gap is between 2 and 
2.5 dB. Two factors contribute to this gap; first, the LDPC 
codes have finite length (coding gap) and second, the input 
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Fig. 1. Comparison of rate-matching by non-equiprobable signaling to 
conventional equiprobable signaling. Both schemes use LDPC codes from 
the DVB-S2 standard. Gains of up to 4 dB (at a spectral efficiency of 
4.6 bits/s/Hz) are obtained. The rate-matched scheme uses two modes, 
namely 64-ASK with a rate 9/10 code and 32-ASK with a rate 5/6 code. 
Equiprobable signaling uses 12 modes, which combine {16, 32, 64}-ASK 
with (3/4,4/5,5/6,8/9,9/10} rate codes. The two schemes are evaluated 
at a frame error rate (FER) of 10“^. For comparison, the operating points of 
the turbo coded modulation schemes with probabilistic shaping suggested in 

(red x) and (red +) are displayed. The operating points of 0.0 
nave a bit error rate (BER) of 10 Since BER < FER, our rate-matched 
scheme is more energy efficient than the schemes suggested in 1^, 0, see 
Sec.liTAl for a further discussion. 

distribution is uniform (shaping gap). At the lower comer 
points, the gap can be as large as 4.5 dB because the system 
resorts to one of the supported transmission modes. 

To provide modes with a finer granularity, we can increase 
the number of code rates and constellation sizes. This approach 
is taken in the extension DVB-S2X Q of DVB-S2. The system 
complexity increases with the number of supported code rates, 
which suggests the use of rate-compatible codes 0. In 0- 
rate-compatible LDPC codes are designed. 

In this work, we take a different approach and propose a 
new coded modulation scheme that uses probabilistic shaping 
to solve the problem of coarse mode granularity and to remove 
the shaping gap. With ASK constellations and only one LDPC 
code rate per constellation size, our scheme operates within 
less than 1 dB of capacity over the whole considered range 
of SNR, see Fig. For the transmitter, we propose a new 
scheme, which we call probabilistic amplitude shaping (PAS). 
At the receiver, bit-metric decoding Gg-© is used. No 
iterative demapping is required. Our scheme directly applies 
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to two-dimensional quadrature amplitude modulation (QAM) 
constellations by mapping two real ASK symbols to one 
complex QAM symbol. 

This work is organized as follows. In Sec. |II] we give an 
overview of the related literature on coded modulation with 


probabilistic shaping. In Sec. Ill we review the information 


theoretic limits of the AWGN channel and we discuss optimal 
signaling for ASK constellations. We introduce PAS in Sec. IV 
and we show in Sec. |V] how PAS can be implemented by 
distribution matching. In Sec. |VI| we combine PAS at the 
transmitter with bit-metric decoding at the receiver. We discuss 


in Sec. VII the code design problem posed by our scheme and 
present a bit-mapper optimization heuristic. We propose a rate 


adaption scheme in Sec. VIII and we discuss numerical results 
in Sec. m The work is concluded in Sec. 


II. Related Literature 

On the additive white Gaussian noise (AWGN) channel, uni¬ 
formly distributed inputs can be up to 1.53 dB less power effi¬ 


cient than Gaussian inputs 113 Sec. IV.B], see also Sec. III-D 


of this work. Coded modulation with probabilistic shaping 
uses a non-uniform distribution on equidistant signal points to 
overcome the shaping gap. For works with a theoretical focus 
on the topic, see for example mm and references therein. 
Here, we review some of the existing works that address the 
practical implementation of probabilistic shaping. 

A. Callager’s Scheme 


Gallager proposes in 116 p. 208] to use a many-to-one 


mapping to turn a channel with non-uniform inputs into a 
super channel with uniform inputs. This approach is combined 
with turbo coding in Q and |j^. In i fTT] , Gallager’s scheme 
is combined with convolutional codes. Optimal mappings are 
investigated in p8) , see also m Sec. LB]. When Gallager’s 
scheme is combined with binary forward error correction 
(FEC), undoing the many-to-one mapping at the receiver is 
challenging. In 0, 0, (Tt), this inversion is achieved by 
iterative demapping, which increases the system complexity 
compared to uniform signaling. In 0, the authors choose the 
many-to-one mapping according to the desired transmission 
rate. The resulting rate granularity is coarse. To achieve a finer 
granularity, more bits must be mapped to the same signal point, 
which increases the system complexity further. 


B. Trellis Shaping 

Trellis shaping is proposed in pO) and it is discussed in 
detail in |21 Chap. 4] and |22 Sec. 4.4]. The transmitter first 
selects a set of sequences. A shaping code then selects from 
this set the sequence of minimum energy for transmission. 
The authors of p3) combine multilevel coding with trellis 
shaping using convolutional codes, see als o p4] Sec VIII.]. 
Bit-interleaved coded modulation (BICM) is used 

in m for lower bit levels and a convolutional shaping code 
selects the two highest bit-levels for energy minimization. 


LDPC codes were used as shaping codes in 1281. Since shaping 
is achieved by the decoder of a shaping code, it is difficult to 
make trellis shaping flexible to support different shaping rates. 


C. Shell Mapping 

Shell mapping | |29) was proposed independently by several 
groups around the same time, see [SO] Sec. VIL], pT| Chap. 8]. 
In shell mapping, the transmitter first selects a shell, which 
is a set of low energy input sequences. It then selects a 
sequence within the selected shell for transmission. Shell 
mapping combined with trellis coding is used in the ITU- 
T Recommendation V.34 m- The main challenge is to 
efficiently index the shells. 


D. Superposition Coding 


Superposition coding was proposed in |32|. The transmitter 


uses several encoders. Each encoder maps its output to a 
signal point. The selected signal points are then superposed 
and transmitted. The receiver uses multistage decoding. Turbo 
codes are considered in and variations of superposition 
coding are proposed and analyzed in p^ . The author of 
0 uses LDPC codes and proposes signal constellations 
that reduce the number of decoding stages. Superposition 
coding requires sequential decoding of the individual codes 
or iterative demapping. 


E. Concatenated Shaping 

The authors of | [35) concatenate BICM encoding with 
shaping. A turbo encoder is followed by a one-to-one shaping 
block code. The scheme is improved in by iterative 
demapping. A variation of this scheme is presented in pT) , 
where the one-to-one shaping block code is followed by 
a many-to-one mapping reminiscent of Gallager’s scheme. 
The shaping scheme p^ is used in p8) with LDPC codes 
on amplitude and phase-shift keying (APSK) constellations. 
Since the shaping decoder is placed before the PEC decoder, 
only limited gains are observed in p5) . The works p^-pS] 
therefore use iterative demapping. 


E Bootstrap Scheme 


A bootstrap scheme is proposed in pg, ig Chap. 7], 
which separates shaping and PEC by first shaping the data and 
then systematically encoding the shaped data. The generated 
check bits are also shaped and embedded in the shaped 
data of the next transmission block. The bootstrap scheme is 
implemented in p9) , [|40| Chap. 7] using Geometric Huffman 


Coding 1411 for shaping and LDPC codes from the DVB- 
S2 standard for PEC. Since shaping is done prior to PEC 
encoding, the bootstrap scheme borrows from the reverse con¬ 
catenation proposed in | |42) , | [43) for magnetic recording. The 
authors in m call the bootstrap scheme chaining construc¬ 
tion and prove that it is capacity-achieving for any discrete 
memoryless channel. The drawback of the bootstrap scheme 
is that shaping is done over several consecutive transmission 
blocks, which have to be decoded in reverse (bootstrapped). 
This increases the latency and is prone to error propagation 
over blocks. 
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TABLE I 

SNR GAPS IN DB of X* TO UNIFORM ASK AND CAPACITY C(P). 


constellation 

rate [bits/channel use] 

X* 

SNR [dB] 

uniform ASK 

SNR [dB] gap [dB] 

capacity C(P) 

SNR [dB] gap [dB] 

4ASK 

1 

4.8180 

5.1181 

-0.3001 

4.7712 

0.0468 

8ASK 

2 

11.8425 

12.6187 

-0.7761 

11.7609 

0.0816 

I6ASK 

3 

18.0910 

19.1681 

-1.0772 

17.9934 

0.0975 

32ASK 

4 

24.1706 

25.4140 

-1.2434 

24.0654 

0.1052 

64ASK 

5 

30.2078 

31.5384 

-1.3307 

30.0988 

0.1090 


III. Optimal Signaling over the AWGN Channel 
The discrete time AWGN channel at time instance i is 
described by the input-output relation 

Y, = X, + Z, (2) 

where the noise terms Zi, i = 1,2,3,... are independent and 
identically distributed (iid) according to a zero mean Gaussian 
distribution with variance one. For i\ channel uses, the input 
is subject to the power constraint 

EEr..|A.I^] ^ p (3, 

ric 

where E[-] denotes expectation. 


C. Optimization of the ASK Input 

We use signal point x d X with probability 


PxAx) = A,= 




( 8 ) 


The scalar ensures that the probabilities assigned by Px,, 
add to 1. Distributions of the form (S are called Maxwell- 
Boltzmann distributions, see, e.g., | |44| . For a fixed constella¬ 
tion scaling A and power constraint P, we choose the input 
distribution 


PxA^) =PxAx) with I.: E[|AA,|2] = P. (9) 


A. AWGN Capacity 

Let X'^ = Xi, X 2 , ■ • ■, be the output of the decoder 
that estimates the input X^ from the output Y"^. The block 
error probability is Pe = PrjA'^ X^). The channel coding 
theorem Theorem 7.4.2] states that by choosing 1 % large 
enough, P^ can be made as small as desired if the transmission 
rate R is smaller than the channel capacity C(P). The goal 
of this work is to design a modulation scheme that achieves 
reliable transmission close to the capacity-power function Q, 
i.e., for any average signal power P, we want to reliably 
transmit data over the AWGN channel at a rate that is close 
to C(P). 

B. Amplitude Shift Keying 

Let X be distributed on some finite alphabet X of signal 
points and suppose E[|Ap] = P. The channel coding theorem 
states that reliable transmission at rate R with average power 
P is possible if 

R<l{X-,Y) (4) 

where !(•; •) denotes the mutual information in bits. The first 
step in designing a coded modulation system is thus to choose 
an alphabet X and a distribution Px on X such that I(A; Y) 
is close to the maximum value C(P). As input alphabet, we 
use an ASK constellation with 2"^ signal points, namely 

T' = {±1,±3,...,±(2”^-1)}. (5) 

We scale X by the constellation scaling A > 0 and the 
resulting input/output relation is 

Y = AX+ Z. (6) 

The power constraint for X is now 

E[|AAp] < P. (7) 


The distribution Px^ maximizes entropy subject to the power 
constraint 0. see | |45| Chapter 12]. Some basic manipulations 
show that E[|Xi,p] is strictly monotonically decreasing in ly. 
Thus, the value of for which the condition (|^ is fulfilled 
can be found efficiently by using the bisection method. For 
each constellation scaling A, the distribution Px^ satisfies the 
power constraint. We now maximize the mutual information 
over all input distributions from this family, i.e., we solve 

maxI(AA;AAA-f ^). (10) 

The mutual information I(Aa;AAa + Z) is a unimodal 
function of A and the optimization problem can be solved 
efficiently by using the golden section method. We denote the 
maximizing scaling by A* and the corresponding distribution 
by Px*. 


D. Shaping Gap 

In Table [1] we display the shaping gains of our input X* 
as compared to a uniformly distributed input. For increasing 
rates and constellation sizes, the shaping gain increases and 
approaches the upper bound of lOlogj^g ^ 


1.53 dB 113 


Sec. IV.B]. The bound can also be derived by probabilistic 
arguments, see | |46l Comment 4)]. For 64-ASK and a rate of 
5 bits per channel use, the shaping gain is 1.33 dB. 


Remark 1. The distribution Px* A suboptimal in general. 
The optimal distribution Px* can be approximated numer¬ 
ically by using the Blahut-Arimoto algorithm 
However, there is no anal^ical expression of Px*. For the 
operating points in Table [d the energy efficiency of A* is 
within 0.1 dB of capacity C(P), so the gain from optimizing 
further is bounded by 0.1 dB. 
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E. Input Distribution: Operational Meaning 


Suppose we use the channel ric times, i.e., we transmit 
length rt codewords. We use codeword with probability 
Px”c If our transmission rate R is larger than the average 
mutual information between input and output, i.e., if 


R > 


i\ 


(11) 


then the probability of error Pr(X" ^ X”) is bounded 
away from zero for any decoder. This result follows as a 


corollary of the converse of the channel coding theorem 116 
Theorem 7.3.1]. 

The mutual information terms l{Xi\Yi) are calculated 
according to the marginal distributions Pxi- We call the 
condition o the Coded Modulation Converse. The insight 
we take from o is that building a transmitter with marginal 
input distributions 


Pxi ^ Px* 


( 12 ) 


Dk 


Dk 


-{P]-> R"->‘ 

Fig. 2. Uniform check bit assumption Sec. 7.1.3]: Data bits = 
D 1 D 2 • ■ ■ Dk with some arbitrary distribution Pok (the data bits are pos¬ 
sibly dependent and non-uniformly distributed) are encoded by a systematic 
generator matrix [J/d-f*] of binary code. The encoder copies the 

data bits to the output, which preserves the probability distribution of 
D^. The encoder calculates n — k redundancy bits by multiplying 

with the parity matrix P. Since each Ri is a modulo two sum of several data 
bits, we assume that the redundancy bits are approximately uniformly 

distributed. 


modulo two sum of a large enough number of bits is close 
to uniform, by a central limit theorem kind of argument. The 
uniform check bit assumption |40 Sec. 7.1.3] states that the 
marginal distributions Pn- are uniform for a large class of data 
distributions Pda:. 


is necessary to achieve reliable transmission close to the rates 

displayed for X* in Table ^ g Encoding Procedure 


IV. Probabilistic Amplitude Shaping 
We now develop a transmitter with property ( [T^ . 

A. Preliminaries 

We make the following two observations: 

1) Amplitude-Sign Factorization: We can write V* as 

X^^A-S (13) 

where A = \X*\ is the amplitude of the input and where 
S = sign(Ar*) is the sign of the input. By Q, the amplitudes 
take values in 

.4:={1,3,...,2™-1}. (14) 

We see from ([^ that the distribution Px* is symmetric around 
zero, i.e., we have 

Px*{x) = Px*i-x) (15) 

and therefore, A and S are stochastically independent and S 
is uniformly distributed, i.e., we have 

= fUda;!) • Ps(sign(a;)), \/x G X (16) 

Ps{l) = Psi-l) = I- (17) 

2 ) Uniform Check Bit Assumption: The second observation 
is on systematic binary encoding. A systematic generator 
matrix of an (n, k) binary code has the form 

G=[h\P] (18) 

where Jj. is the k x k identity matrix and P is a fc x (n — 
k) matrix. P is the parity matrix of the code. The generator 
matrix G maps k data bits to a length n codeword via 

D'^G = (D^R^-^) (19) 

where are redundant bits that are modulo-two sums of 

data bits. See Fig. for an illustration. The distribution of the 


Consider block transmission with symbols from a 2"*- 
ASK constellation. Since we use binary error correcting codes, 
we label each of the 2"*“^ amplitudes by a binary string of 
length TO — 1 and we label each of the signs ±1 by a bit, i.e., 
we use 


Ah^6(A) e {0,1}™-^ (20) 

S ^biS) G{0,1}- ( 21 ) 


For the sign, we use h(—1) = 0 and 5(1) = 1. We discuss the 
choice of b{A) in Sec. VTC We use a rate k/n = {m — l)/m 
binary code with systematic generator matrix G — [7fc|P]. 
For block transmission with 1 % channel uses, the block length 
of the code is n = rfcTO and the dimension of the code is 
k = The encoding procedure is displayed in Fig. 

It works as follows. 


1) A discrete memoryless source (DMS) \Pa\ outputs 
amplitudes Ai, A 2 ,..., A^, that are iid according 
to Pa- We explain in Sec. [v] how the DMS \Pa\ 
can be emulated from binary data by distribution 
matching. 

2) Each amplitude Ai is represented by its label 

b{Af. 

3) The resulting length (to — 1)7^ = k binary string 
is multiplied by the parity matrix P to generate 
n - k = nc sign labels 5(S'i), 5(52), ■ ■ ■, 6(5,^). 

4) Each sign label b{Si) is transformed into the 
corresponding sign Si. 

5) The signal Xi = Ai ■ Si is scaled by A and 
transmitted. 

We call this procedure probabilistic amplitude shaping (PAS). 
Since the signs S"^ are a deterministic function of the ampli¬ 
tudes A^, the input symbols Vi, ^ 2 ,..., X^^. are correlated. 
Let’s check if the marginal distributions Px^ fulhll condition 
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Fig. 3. PAS. The ASK amplitudes Ai take values in ^ = {1,3,..., 2^ — 1}. The amplitudes Ai ars represented by their binary labels b{Ai). Redundancy 
bits b{Si) result from multiplying the binary string b(Ai)b{A 2 ) • • ■ b{Ari,) by the parity matrix P of a systematic generator matrix [J|P]. The redundancy bits 
b{Si) are transformed into signs Si and multiplied with the amplitudes Ai. The resulting signal points Xi = AiSi take values in A = (±1, ±3,..., ±(2^^ — 
1)}. The signal points Xi are scaled by A and AXi is transmitted over the channel. 


0- We have 


Pxiixi) = P^,Si(|xi|,sign(a;i)) (22) 

^ PAi{\xi\)Psiisign{xi)) (23) 

= PA(|a;*|)Psi(sign(xi)) (24) 

^Pa{\x^\)1 (25) 

= Px*ix) (26) 


where we used the uniform check bit assumption in ( [25] l. We 
discuss when we have equality in ( [23] l and ( [25] l. 

• We have equality in ( |2^ if Ai and Si are independent. 
This is the case if the redundant bit b(Si) does not depend 
on the binary label of Ai (the bit b{Si) can still depend 
on other Aj, j ^ i). On the other hand, if b{Si) is 
determined by the binary index of Ai only, then Ai and 
Si are dependent and we will not have equality in ( [23] l. 

• Recall that h{Si) is the modulo two sum of binary random 
variables. If one of these binary random variables is 
independent of the others and uniformly distributed, then 
a basic probability calculation shows that b{Si) is exactly 
uniformly distributed and we have equality in ( |25] l. 

For the rest of this work, we assume that our scheme achieves 
equality in ( |2^ and ( |25| ). Our numerical results support this 
assumption. Of course, if the empirical results deviate from 
what theory predicts, we should check the assumptions of 
equality in (|23|l and (|25]l. 



SNR in dB 


Fig. 4. The mutual information curves (solid) and the transmission rate curves 
(dashed) for ASK. The optimal operating points for rate (m — l)/m codes 
are indicated by dots. 

C. Optimal Operating Points 

We study the rates at which reliable transmission is possible 
with our scheme. By 0, reliable communication at rate R 
is achievable only if 

R < = i(X; Y) = I(y45'; Y). (28) 


Remark 2. PAS is a special case of the bootstrap scheme 
0 0- After encoding, the redundancy bits are uniformly 
distributed. However, instead of transforming them into a 
sequence of symbols with a non-uniform distribution and 
transmitting them in the next block as is done in the bootstrap 
scheme, they can be used directly as sign labels in the same 
block, since the uniform distribution is already optimal. 

Remark 3. PAS can be seen as a probabilistic version of shell 
mapping. By the weak law of large numbers, 

rh a2 

^ « P (27) 

tie 

with high probability, i.e., the source Pa selects an amplitude 
sequence in a shell in the rtc-dimensional space that contains 
the sequences of energy « rteP. Multiplying the amplitudes by 
signs does not change the power, i.e., the transmitted signal 
remains in the shell selected by the amplitude source Pa- 


Since A^ represents our data, our transmission rate is 


ric 


H(A) 


and condition ([28]) becomes 


bits 

channel use 


(29) 


H(y4) < 1(^5;^). (30) 

In Fig. 1^ both the mutual information I(AS';y) (solid lines) 
and transmission rate IHI(A) (dashed lines) are displayed for 
4, 8,16, 32, and 64-ASK. For high enough SNR, the mutual 
information saturates at m bits and the transmission rate 
saturates at m — 1 bits. Optimal error correction for block 
length Uc —oo would operate where the transmission rate 
curve crosses the mutual information curve. These crossing 
points are indicated by dots in Fig. How to operate at other 
transmission rates is the topic of rate adaption and we discuss 
this in detail in Sec. Eul 
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Fig. 6 . Extension of PAS to code rates higher than (m— l)/m. The fraction 7 of the signs is used for data, which is modelled as the output of a Bemoulli-1/2 
DMS Pu- 



Fig. 5. Optimal operating points of 8 -ASK for PAS (c = 2/3) and extended 
PAS (c = 3/4). 


D. PAS for Higher Code Rates 

We observe in Fig. that the ASK mutual information 
curves stay close to the capacity C(P) over a certain range of 
rates above the optimal operating points. We therefore extend 
our PAS scheme to enable the use of code rates higher than 
(to — 1)/to on 2”^-ASK constellations. We achieve this by 
using some of the signs Si for uniformly distributed data bits. 
We illustrate this extension of the PAS scheme in Fig. Let 7 
denote the fraction of signs used for data bits. We interpret 
uniformly distributed data bits as sign labels b{Si) ■ ■ ■ 6(S'-y^). 
These 777 bits and the (to — 1)77 bits from the amplitude 
labels are encoded by the parity matrix of a rate c code, which 
generates the remaining (1 — 7 ) 7 ^ sign labels. The code rate 
can be expressed in terms of to and 7 as 


TO — 1 + 7 

C = - 

TO 


(31) 


For a given code rate c, the fraction 7 is given by 


7 = 1 — (1 — c)to. 


(32) 


Since a fraction 7 of the signs now carries information, the 
transmission rate of the extended PAS scheme is given by 


rh 


bits 


channel use 


(33) 


UiU2...Uf^ 


Matcher 


A1A2... Ar^ 


Fig. 7. The matcher transforms uniform data blocks of length At into 
rt amplitudes that are approximately distributed according to the desired 
distribution By replacing the amplitude source in the PAS diagrams 
in Fig. and Fig. ||]by a matcher, our scheme provides a binaiy interface to 
the source coding part of a digital communication system. 


The optimal operating point is then given by the crossing of 
the rate curve IHI(A) +7 and the mutual information. In Fig. 
we display for 8 -ASK the optimal operating points for c = 2/3 
and c = 3/4. 


V. Distribution Matching 


In a digital communication system, usually a binary inter¬ 
face separates the source coding part from the channel coding 
part [ |49l Chap. 1]. Up to now, our scheme does not have such 
a binary interface, since some of our data is output by an 
amplitude source \Pa\ - We therefore add a device that takes as 
input uniformly distributed independent bits and outputs 
77 c amplitudes A’^. We illustrate this in Fig. We require the 
following properties for our device: 

PI The input provides a binary interface to the source coding 
part of the digital communication system. 

P2 The output is approximately the output of a DMS IP 4 I . 

P3 The input can be recovered from the output, i.e., the 
mapping is invertible. 


A device with such properties is called a distribution matcher 


|41|. Variable length distribution matchers were proposed in 
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Sec. IV.A], ||44| Sec. VII] and fSO) and their design 


studied in |g, 

Variable length matchers can lead to buffer overflow, syn¬ 
chronization loss and error propagation, see, e.g., |44 Sec. 1]. 
We therefore use the fixed length distribution matcher pro¬ 
posed in 0 . This matcher, called constant composition 
distribution matcher (CCDM), has the properties PI and P3. 
To address property P2, we need to say what we mean by 
“matcher output A^ is approximately the output of a DMS 
I Pa I ”. By Theorem 1.2], a good measure for similarity in 
the context of channel coding is the normalized informational 
divergence (also known as Kullback-Leibler divergence or 
relative entropy | |45| Sec. 2.3]) of the output distribution 
and the memory less distribution P^- The normalized 
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informational divergence is 




E 

Gsupp 




P^^(a )10g2 


rh 


where suppP^^ is the support of P^,^. The CCDM has 
property P2 in the following sense: 

1) As rfc —> oo, the normalized informational divergence 
( |T4| i approaches zero. 

2) As r^; —>■ oo, the rate approaches IHI(A), i.e., we have 

-^H(A). (35) 

rh 

We use the CCDM to emulate the amplitude source m The 
output length is finite, and we use the actual rate for 
performance evaluations. 


VI. Bit-Metric Decoding 
The receiver estimates the transmitted codeword from 
the channel outputs . In this section, we show how this can 
be implemented by a bit-metric decoder. 


A. Preliminaries: Binary Labeling and Decoding 

Consider a coded modulation system where the input X 
takes values in a 2”^-ASK constellation. An optimal decoder 
uses the symbol-metric py\x 111 Sec. II.B] and can achieve 
the rate 


RsMD=I(^;n (36) 

where SMD stands for symbol-metric decoding (SMD). We 
are interested in successfully decoding at a transmission rate 
close to Rsmd by using a binary decoder. Recall that our PAS 
scheme labels the amplitude A = | ATj by a length to — 1 binary 
string b{A) and the sign S = sign(2f) by one bit b{S). The 
length TO binary string 

B = BiB2---B^:=b{S)biA) (37) 

thus assigns to each signal point x G X a label via 

label(a:) = label(sign(a:))label(|a;|) = 61&2 • • • b^. (38) 

Since the labeling is one-to-one, we can also select a signal 
point for transmission by choosing the label B, i.e., we use 

X = XB = {x G X:lahA{x) = B}. (39) 

We can interpret B as the channel input and the input/output 
relation of our channel as 


Y = Axb + Z. 


(40) 


Using the chain rule, we expand the mutual information of 
B and Y as 


= I{ByY)+I{B2-,Y\Bi) + --- 

+ Y\Bi ■ ■ ■ Bm-l)- 


(41) 


This expansion suggests the following binary decoding: 

1) Use the channel output Y to calculate an estimate Bi. 

2) Successively use the output Y and the estimates 
Bi ■ ■ ■ Bi-i to calculate the next estimate Bi. 

This approach is called multistage decoding (MD). It requires 
multilevel coding (MLC) at the transmitter, i.e., on each bit- 
level, an individual binary code with block length rfc is used. 
MLC/MD was first introduced in ||5^ and it is discussed in 


detail, e.g., in |24|. To use MLC/MD, we would need to 
modify our PAS scheme. A simpler approach is to ignore the 
estimates Bj, j ^ i when estimating Bi. This reduces the 
mutual information, which can be seen as follows 


m 

I{B;Y)=Y,HBy,Y\B^-^) 

m 

= '^E.{Bi\B^-^) -E.{Bi\YB^-^) 

(a) ^ 

m 

= H(B) 

2=1 


(43) 

(44) 

(45) 

(46) 


where (a) follows because conditioning does not increase 
entropy |45 Theorem 2.6.5] and where we used the chain 
rule for entropy in (b). The expression in the last line of 
(gg can be approached as follows. We jointly encode all 
bit-levels by a single binary code of block length rm\. This 
idea was introduced in | [25) and is now usually called bit- 
interleaved coded modulation (BICM) Since our 

PAS transmitter encodes all bit-levels by a single binary code 
G = [/fc|P], our transmitter is a BICM encoder. At the 
receiver, we use a bit-metric decoder. 


B. Bit-Metric Decoding 

A soft-demapper calculates for each bit-level i the soft- 
information 


T^ = 



+ log 


PY\Bm^) 

PY\Bi{X\Yj 

- ^ 


a-priori information channel likelihood 


(47) 


The distribution and the conditional density PY\Bi in 
can be calculated as 


PB,{bi) = ^ PB{a) 

aG{0,l}^: ai—bi 

PY\BM^i)= E PY\B{v\a) 


(48) 

(49) 


The soft-information P, forms a sufficient statistic to estimate 


bit level Bi from the channel output Y |57 Sec. 8.2.1], i.e., 
we have I{Bi;Y) = 1(5^; Li). A bit-metric decoder uses the 
soft information Li, L 2 ,..., Lm to estimate the transmitted 
data and achieves the rate ifTT] Theorem 1] 


Rbmd — H(P)-^H(P,|r). (50) 

i=l 


( 42 ) 













TABLE II TABLE III 

Two LABELINGS FOR THE AMPLITUDES OF 8-ASK AND THE RESULTING COMPARISON OF SNRS NEEDED BY SYMBOL-METRIC DECODING AND 

SIGNAL POINT LABELINGS BIT-METRIC DECODING TO ACHIEVE A CERTAIN RATE 




Amplitude 

7 

5 

3 

1 



constellation 

rate 

SNR(Rsmd) [dB] 

SNR(Rbmd) [dB] 

Gap [dB] 


natural labeling 

00 

01 

10 

11 



4-ASK 

1 

4.8180 

4.8313 

0.0133 



BRGC 

00 

01 

11 

10 



8-ASK 

2 

11.8425 

11.8481 

0.0056 

Signal point 

-7 

-5 -; 


-1 

1 

3 

5 

7 

16-ASK 

3 

18.0911 

18.0951 

0.0039 

Labeling 1 

000 

001 010 

oil 

111 

no 

101 

100 

32-ASK 

4 

24.1708 

24.1742 

0.0034 

BRGC 

000 

001 oil 

010 

110 

111 

101 

100 

64-ASK 

5 

30.2078 

30.2110 

0.0032 


SNR in dB 



Eig. 8. Comparison of Rbmd for the two labelings from Table [Ti] 

Note that ( |50l l is equal to the expression we derived in ( |46l l. 

Remark 4. If the bit levels are independent, the rate becomes 
the “BICM capacity” €3 Theorem 1 ] 

m 

(51) 

i=l 

Using non-uniform distributions for the independent bit-levels 
results in bit-shaping, which was introduced in H58^ and is 
discussed in & The authors in / [59) Eq. (1)] call (HB the 
“pragmatic capacity”. 


is a natural labeling and Labeling 2 is the Binary Reflected 
Gray Code (BRGC), see | |60) . The two labelings are listed 
in Table |I^ We have ordered the amplitudes in descending 
order so that the actual signal points are labeled from left 
to right. The resulting labelings of the signal points are also 
displayed in Table |II] Note that the natural labeling of the 
amplitudes does not lead to a natural labeling of the signal 
points. In contrast, the BRGC labeling of the amplitudes also 
leads to the BRGC of the signal points. We display in Fig. 
Rbmd for Labeling 1 and the BRGC labeling. The BRGC 
labeling is better than Labeling 1 and very close to the capacity 
C(P), which is consistent with the results presented in p^. 


In Table III bit-metric decoding with BRGC is compared to 
symbol-metric decoding. The SNRs needed to achieve a given 
rate are listed. 


Remark 5. The bit-metric decoding gaps in Table show 
how much we can gain by using multilevel coding with 
multistage decoding or by iteratively exchanging extrinsic 
information between the soft-demapper and the soft-decoder 
(BICM-ID, H61]l j. This gain is negligible and not worth the 
increased complexity for the considered scenario. 


VII. LDPC Code Design 

In principle, our scheme works for any binary code with a 
systematic encoder and a decoder that can process the soft- 
output of the binary demappers. In this section, we discuss the 
deployment of LDPC codes. 


A. LDPC Codes and Bit-Channels 


From a coding perspective, bit-metric decoding transforms 
the channel ( |40| l into m parallel binary input channels 
i = 1, 2,..., m. 


are different for different bit 


1) The bit channels PLi\Bi 
levels. This is important in Sec. |VII| where we optimize 
the parity matrix P for bit-metric decoding. 

2) By ( |49l l, the distribution of the bits Bj, j i influences 
the channel transition probabilities PLi\Bi of the bit- 
channel, i.e., for each bit level, the other bit-levels act 
as interference. 

3) The rate (|^ can take on different values for different 
labelings of the signal points. We discuss this next. 


C. Optimizing the Labeling 

We choose the labeling of the to — 1 amplitudes; the label of 
the sign is already defined by construction. We evaluate Rbmd 
for 8-ASK and two different amplitude labelings. Labeling 1 


LDPC codes are linear block codes with a sparse {n—k)xn 
check matrix H. The matrix H can be represented by a Tanner 
graph | [62| consisting of variable nodes and check nodes. The 
variable node degree of the ith coded bit is given by the 
number of ones in the Ah column of H and the degree of 
the jth check node is given by the number of ones in the jth 
row of H. The variable and check node degrees influence the 
decoding threshold of the LDPC code ||6^. 

Good LDPC codes are often irregular, i.e., not all coded bits 
have the same variable node degree | [63] , | [64| . At the same 
time, the coded bits are transmitted over different bit-channels. 
This suggests that the bit-mapper, which decides which coded 
bit is transmitted over which bit-channel, influences the per¬ 
formance. The example shown in Fig. confirms this. The 
optimized bit-mapper is 0.6 dB more energy efficient than a 
randomly chosen bit-mapper. 

Bit-mapper optimization was considered, e.g., in |[^-|[^. 
An alternative approach is to jointly optimize the node degrees 
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SNR in dB 

Fig. 9. Influence of the bit-mapping on the performance. The DVB-S2 rate 
2/3 LDPC code is used with uniform inputs on an 8-ASK constellation. An 
optimized bit-mapper (blue curve) and a random bit-mapper (red curve) are 
used. 


and the bit-mapper. This is done in fTO) and | |7Tj . In this work, 
we consider bit-mapper optimization. 


B. Bit-Mapper Optimization for the DVB-S2 Codes 

We use LDPC codes from the DVB-S2 standard. In Ta¬ 


ble IV we display the variable node degree distributions of 
the DVB-S2 codes, e.g., the rate 2/3 code has 4320 coded bits 
of variable node degree 13, and it has 38880, 21599, and 1 
bits of degrees 3, 2, and 1, respectively. All codes have four 
different variable node degrees, which decrease from left to 
right. By Sec. IV-B our PAS scheme places the uniformly 
distributed bits of bit-level Bi at the end of the codeword in 
the systematic encoding process. For the bit-levels B 2 ■ ■ ■ Bm, 
we use the following heuristic. 

• A bit interleaver TTb sorts the bit stream by bit-levels, i.e.. 


6(Ai)...6(A^) 


BrfBy. 


B„ 


(52) 


where Bi is a string of i\ bits of level i. 

• A bit-level interleaver permutes the bit-level strings. We 
compactly represent the bit-level interleaver by listing 
the bit-levels in the order in which they should occur 
in the codeword, e.g., for m = 4, the bit-level interleaver 
(4,2,3) is 

B 2 -B 3 B 4 B 4 B 2 S 3 . (53) 


We keep the bit interleaver fixed and represent the complete 
bit-mapping by appending a 1 to the bit-level interleaver, e.g., 
TT = (4, 2, 3,1) stands for the composition of TTb with bit- 
level interleaver (4,2, 3). We optimize the bit-level interleaver. 
There are {m — 1)! possibilities, among which we choose the 
one with the best error performance. The largest considered 
constellation is 64-ASK, for which we need to choose among 
(6—1)! = 120 bit level interleavers. This optimization is still 
feasible. We display in Table [V] the optimized bit-mappers for 
uniform inputs and rate (m — l)/m codes. 


TABLE IV 

Variable Node Degree Distributions of DVB-S2 Codes 


rate 



vaiiable node degrees 



13 

12 

II 

4 

3 

2 

1 

2/3 

4320 




38880 

21599 

1 

3/4 


5400 



43200 

16199 

1 

4/5 



6480 


45360 

12959 

1 

5/6 

5400 




48600 

10799 

1 

9/10 




6480 

51840 

6479 

1 


TABLE V 

Optimal bit-level interleaver for uniform inputs and DVB-S2 

CODES 


constellation 

code rate 

bit-mapper tt 

8ASK 

2/3 

(3,2,1) 

16ASK 

3/4 

(3,4, 2,1) 

32ASK 

4/5 

(3, 2, 5, 4,1) 

64ASK 

5/6 

(4, 2, 5, 3,6,1) 


VIII. Rate Adaption 
A. Practical Operating Points 

We use PAS at the transmitter and bit-metric decoding at 
the receiver. The transmission rate and the achievable rate are 
given by the respective 


i? = H(A)+ 7 , RBMD=H(B)-^iT(B,|L0 (54) 

i=l 


where B = b{S)b{A). For each transmit power P, we choose 
the amplitude distribution Pa and the constellation scaling 
A to maximize the achievable rate. In this way, we obtain 
a transmission rate curve R and an achievable rate curve 
Rbmd- In analogy to our discussion in Sec. IV-C the optimal 
operating point for our scheme is where R crosses Rbmd- We 
illustrate this in Fig. 10 for 8 -ASK and code rate c = 2(3 
(7 = 0). When we use practical codes of finite block length 
Tic, we must back off from the optimal operating point and 
tolerate a positive FER. We back off along the transmission 
rate curve R. As we increase the transmission power P, the 
rate back-off 


R 


BMD 


- R 


(55) 


increases and we expect that the error probability decreases. 
This intuition is confirmed by the practical operating points of 


the 2/3 DVB-S2 LDPC code shown in Fig. 10 


B. Transmission Rate Adaption 

Suppose we achieve a transmission rate R° at a FER Pe 
by the procedure described in Sec. VIII-A Suppose further 
that this reference operating point is achieved by amplitude 
distribution Pa and constellation scaling A. With the same 
code, we now want to transmit at a rate Rf^R° and achieve 
the same error probability. 
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Fig. 10. Operating points for 8-ASK and rate 2/3 codes. The optimal point 
is where the rate curve R{P) crosses the achievable rate curve Rbmd(P)- 
By backing off from the optimal point along the rate curve, the rate back-off 
(given by the distance between the rate curve and the achievable rate curve) 
increases. For the rate 2/3 DVB-S2 LDPC code, the FER Pg is decreased 
from 1.1 X 10“^ to 2.8 x 10“^ by increasing the rate back-off. 


Recall that for some optimized ly (see Sec. |lII-C] l, the input 
distribution Px* is given by 

Px*ix) = ^ - , xGX. (56) 

By we can write the amplitude distribution Pa of our 
reference point as 

-Pa(o) = 2-Px*(a), a G Jl. (57) 


Now dehne 


Pa^ (a) 


PA(a)e^“" 


a G .A. 


(58) 


The distribution has the following properties; 

• Pa and are Maxwell-Boltzmann distributions. In 
particular, the set of Maxwell-Boltzmann distributions is 
closed under the mapping Pa i—> Pa^ as defined by ( [58] l. 

• For A = 0, we have P^a = Pa- 

• For X ^ V, Pa approaches the uniform distribution on 
A whose entropy is log 2 |Al| = m — 1. 

• For A —>^ —oo, Pa^ approaches the distribution that 
chooses the smallest amplitude with probability 1. The 
resulting entropy is zero. 

From these properties, we see that we can use P^a to adapt 
the transmission rate. The range of feasible rates is 


7 < + 7 <m —1 + 7 . (59) 


For 7 = 0, the rate is between 0 and m — 1. To transmit at a 
feasible rate R, we proceed as follows. 


• Choose the amplitude distribution such that the 
desired transmission rate is achieved, i.e., choose 

Px = Pa^ ■■ + 7 = P. (60) 


* Choose the constellation scaling A such that the 
resulting error probability Pg is the same as for the 
reference point, i.e., choose 

Pe = Pe- (61) 


C. Adaption for Universal Codes 

The procedure described in the previous section requires 
to search for the right constellation scaling A by repeatedly 
performing Monte Carlo simulations. We now show how this 
can be simplified significantly for universal codes. We say a 
code is universal, if Pe depends only on the rate back-off 
Rbmd ~ R- Suppose for Pa and A of our reference point, 
the achievable rate evaluates to Rbmd- We can now adapt the 
constellation scaling as follows. 

* Choose the constellation scaling A such that 

Rbmd (A, P^) - P = Rbmd " R° ■ (62) 

Remark 6. For PAS and bit-metric decoding, we can write 
the rate back-off as 

ra 

RBMY,-R=l-l-Y,H{B,\Li). (63) 

i=l 

The term 1 — 7 /s constant, so Pe is determined by the sum 
of the conditional entropies R{Ri\Ri)- 

Remark 7. For a uniformly distributed binary input B and a 
rate c code, the transmission rate is R = c and the achievable 
rate is R* = I(P; V). The rate back-off becomes 

R* -R = I(P; Y) - c. (64) 

Since c is constant, this shows that Pe is determined by the 
mutual information I(Jf; V). This property was observed for 
practical LDPC codes in 4^- 


IX. Numerical Results 


We assess the performance of our scheme by Monte Carlo 
simulation. For each block, uniformly distributed data bits 
are transmitted in channel uses, yric of the data bits 
[/k _ U 1 U 2 ■ • • are used for bit-level 1 (see Fig. |^. The 
remaining At — 77 ^ data bits are transformed by a CCDM 
matcher (see Sec. 0 into a sequence of amplitudes. 
Encoding is then done according to Fig. 1^ At the receiver, bit- 
metric decoding is performed (see Sec. |VI| i and the amplitude 
estimates are transformed back into At — 77 ^ data bit estimates 
by a CCDM dematcher, see 0. The dematcher output 
together with the other yrii. data bit estimates form the data 
estimate C/k = UiU 2 ‘--Uf^. We estimate the end-to-end 
FER Pr{t 7 k ^ f/k } por each EER estimate, we also provide 
the corresponding 95% confidence interval (Cl). The spectral 
efficiency is P = We use the DVB-S2 LDPC codes and we 
decode with 100 iterations. We optimize the bit-mappers using 


the heuristic proposed in Sec. VII-B Eor each considered 
(constellation, code rate) mode, we hrst determine the practical 
operating point with EER « 1 x 10“^ following the procedure 
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TABLE VI 

64-ASK, 9/10 DVB-S2, tt = (4, 2, 5, 3, 6,1) 


Rate 

SNR [dB] 

Gap [dB] 

EER 

95% Cl 

5.09 

31.80 

1.15 

4.1 ■ 

10-3 

±2.1 ■ 

10-3 

4.98 

31.09 

1.09 

8.0- 

10-3 

±4.0- 

10-3 

4.89 

30.46 

1.05 

6.8- 

10-3 

±3.4- 

10-3 

4.79 

29.84 

1.03 

7.3- 

10-3 

±3.7- 

10-3 

4.69 

29.23 

1.02 

6.4- 

10-3 

±3.2 ■ 

10-3 

4.59 

28.62 

1.02 

3.8 ■ 

10-3 

±1.9- 

10-3 

4.49 

28.01 

1.01 

6.2 ■ 

10-3 

±3.1 ■ 

10-3 

4.39 

27.41 

1.01 

7.6 ■ 

CO 

1 

O 

1—1 

±3.8 ■ 

10-3 

4.29 

26.80 

1.00 

9.6 ■ 

10-3 

±4.8 ■ 

10-3 

4.19 

26.21 

1.00 

1.0- 

10-3 

±5.1 ■ 

10-3 

4.09 

25.61 

1.00 

1.0- 

10-3 

±5.1 ■ 

10-3 

3.99 

25.00 

0.99 

1.0- 

10-3 

±5.0 ■ 

10-3 



TABLE VII 





32-ASK, 5/6 DVB-S2, 

TT = 

(4, 5, 2, 3,1) 


Rate 

SNR [dB] 

Gap [dB] 

EER 

95% Cl 

3.62 

22.60 

0.82 

2.1 ■ 

10-3 

±1.1 ■ 

10-3 

3.99 

25.09 

1.07 

3.3 ■ 

10-3 

±1.7- 

10-3 

3.89 

24.36 

0.94 

1.8- 

10-3 

±9.0 ■ 

10-4 

3.79 

23.69 

0.87 

2.8 ■ 

CO 

1 

o 

1—1 

±1.4- 

10-3 

3.69 

23.05 

0.84 

1.5 ■ 

10-3 

±8.0 ■ 

10-4 

3.59 

22.42 

0.82 

1.6 ■ 

10-3 

±8.0 ■ 

10-4 

3.49 

21.80 

0.81 

1.5 ■ 

10-3 

±8.0 ■ 

10-4 

3.39 

21.19 

0.80 

1.3 ■ 

10-3 

±6.0 ■ 

10-4 

3.29 

20.58 

0.80 

8.0- 

10-"^ 

±4.0- 

10-4 

3.19 

19.98 

0.80 

7.0- 

10-4 

±4.0 ■ 

10-4 

3.10 

19.37 

0.80 

5.0- 

10-4 

±3.0 ■ 

10-4 

3.00 

18.76 

0.80 

7.0- 

10-4 

±4.0 ■ 

10-4 



TABLE VIII 





16-ASK, 5/6 DVB-S2, tt = 

(4,3,2 

1) 


Rate 

SNR [dB] 

Gap [dB] 

EER 

95% Cl 

2.96 

18.40 

0.67 

2.0- 

10-3 

±9.8 ■ 

10-3 

3.00 

18.66 

0.69 

3.0- 

10-3 

±1.4- 

10-3 

2.90 

18.01 

0.65 

2.3 ■ 

10-3 

±1.2- 

10-3 

2.80 

17.38 

0.63 

2.0- 

10-3 

±9.8 ■ 

10-3 

2.70 

16.75 

0.62 

1.6 ■ 

10-3 

±7.8 ■ 

10-3 

2.60 

16.13 

0.61 

1.3 ■ 

10-3 

±6.7- 

10-3 

2.50 

15.51 

0.61 

7.2- 

10-3 

±3.6 ■ 

10-3 

2.40 

14.89 

0.61 

3.5 ■ 

10-3 

±1.8- 

10-3 

2.30 

14.27 

0.62 

2.5 ■ 

10-3 

±1.3- 

10-3 

2.20 

13.64 

0.62 

1.9 ■ 

10-3 

±9.0 ■ 

10-4 

2.10 

13.01 

0.62 

1.3 ■ 

10-3 

±7.0 ■ 

10-4 

2.00 

12.37 

0.63 

5.0- 

10-4 

±2.0 ■ 

10-4 



TABLE IX 





8-ASK, 3/4 DVB-S2, tt 

= (3,2,1) 


Rate 

SNR [dB] 

Gap [dB] 

EER 

95% Cl 

1.85 

11.45 

0.63 

1.5 ■ 

10-3 

±8.0 ■ 

10-4 

2.00 

12.44 

0.69 

3.8- 

10-3 

±1.9- 

10-3 

1.90 

11.75 

0.64 

2.1 ■ 

CO 

1 

o 

1—1 

±1.1 ■ 

10-3 

1.80 

11.08 

0.62 

1.2 ■ 

10-3 

±6.0 ■ 

10-4 

1.70 

10.41 

0.62 

2.3 ■ 

10-3 

±1.2- 

10-3 

1.60 

9.75 

0.62 

1.2 ■ 

10-3 

±6.0 ■ 

10-4 

1.50 

9.07 

0.63 

1.5 ■ 

10-3 

±7.0 ■ 

10-4 

1.40 

8.39 

0.64 

1.2 ■ 

10-3 

±6.0 ■ 

10-4 

1.30 

7.69 

0.65 

2.2 ■ 

10-3 

±1.1 ■ 

10-3 



TABLE X 





4-ASK, 2/3 DVB-S2, tt 

= (2,1) 



Rate 

SNR [dB] 

Gap [dB] 

EER 

95% Cl 

1.13 

6.70 

0.90 

5.2 ■ 

10-3 

±2.6 ■ 

10-3 

1.20 

7.26 

0.94 

2.6 ■ 

10-3 

±1.3- 

10-3 

1.10 

6.45 

0.90 

6.4- 

10-3 

±3.3 ■ 

10-3 

1.00 

5.66 

0.90 

1.4- 

10-3 

±6.8 ■ 

10-3 



Fig. 11. Transmission rate adaption for 8-ASK. The rate 3/4 DVB-S2 code is 
used. The triangles maiL the operating points calculated under the assumption 
that the code is universal. The refe renc e operating point is the triangle on 
the rate IHI(A) + ^ curve. In Table [dc| we display the FER of the adapted 
operating points. For comparison, vteaisplay an operating point of the 2/3 
DVB-S2 code. 


described in Sec. |VIII^ We then use the rate adaption ( |62l l to 
operate the same (constellation, code rate) mode over a range 
of SNRs and spectral efficiencies. The results are displayed 
in Table Vj^ In the first line of each table, the practical 
operating point is displayed and separated from the rest of the 
displayed values by a horizontal line. For 8-ASK, we display 


the resulting operating points in Fig. 11 


1) Spectral Efficiency: In Table |VH—|X] we display the SNR 
gap of our scheme to capacity C(P) =Tlog 2 (l + P/1)- Over 
a range of 1 bit/s/Hz to 5 bit/s/Hz, our scheme operates within 
1 dB of capacity. The rate 5/6 code with 16-ASK has a gap of 
only 0.61 to 0.69 dB, while the rate 9/10 code with 64-ASK 
has a gap of 1 dB. This indicates that the spectral efficiency 
can be improved further by using optimized codes. The work 
IZ3 confirms this. 


2 ) Universality: The results in Table IVIj—|X| show that the 
DVB-S2 codes are universal in the sense of Sec. |VIII-C| 
since over the whole considered range, the resulting FER is 
within the waterfall region. Some configurations show a “more 
universal” behavior than others, e.g., the FER of the rate 3/4 
code with 8-ASK (Table [I^ is almost the same for the spectral 
efficiencies from 1.3 to 2.0 bits/s/Hz, while for the rate 5/6 
code with 16-ASK, the EER changes by almost two orders of 
magnitude. Note that the listed operating points were found 
by ( |6^ . The operating points that result in a too low (too 
high) EER could be corrected by decreasing (increasing) the 
transmission power. 

3) Rate: Since all considered codes have the same block 
length of 64800 bits, the number of channel uses i\ gets 
smaller for larger constellations. As discussed in Sec. [Vj the 
CCDM matcher that we use to emulate the amplitude source 
Pa has a rate close to EI(A) when the output length is large. 
Eor the smaller constellations 4, 8 and 16-ASK, the effective 
rate is equal to IHI(A) with a precision of two decimal places. 
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For the larger constellations 32 and 64-ASK, we observe a 
slight decrease of the effective rate, e.g., instead of the target 
rate 5.00 we observe an effective rate R = 4.98. See also flU 
for a discussion of this phenomenon. We calculate the SNR 
gap with respect to the effective rate. 

4) Error Floor: The DVB-S2 codes are designed to get the 
FER down to a reasonable value and an outer code then lowers 
the error probability further Q. Since the CCDM dematcher 
performs a non-linear transformation, one corrupted symbol 
at the dematcher input can lead to various corrupted bit-errors 
at the dematcher output, so in-block error propagation can 
occur. In our scheme, an outer encoder should therefore be 
placed between the matcher and the inner encoder and the 
outer decoder should be placed between the inner decoder and 
the dematcher. To preserve the amplitude distribution imposed 
by the matcher, a systematic outer encoder has to be used. 

X. Conclusions 

We proposed a practical rate-matched coded modulation 
scheme that adapts its transmission rate to the SNR. At 
the transmitter, a distribution matcher is concatenated with 
a systematic encoder of a binary LDPC code. At the re¬ 
ceiver, bit-metric decoding is used. No iterative demapping 
is required. The reported numerical results show that on the 
complex baseband AWGN channel, any spectral efficiency 
between 2-10 bit/s/Hz can be achieved within 1 dB of capacity 
log 2 (1-fP/l) by using the off-the-shelf DVB-S2 LDPC codes 
of rate 2/3, 3/4, 4/5 and 9/10 together with QAM constellations 
with up to 4096 signal points. Future work should investigate 
rate-matched coded modulation for short block lengths and for 
multiple-input multiple-output channels. 
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